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Abstract. We list some explicit calculations related to a theorem of Clausen 
originally published in 1828, more commonly known as the result that describes 
the linear third order differential equation satisfied by the squares and the 
product of any two solutions of a linear second order differential equation in 
the real domain. The case of the cube of a solution dates to Appell, 1880. 
Although not commonly known and perhaps even new to some extent we 
show that the fourth and fifth powers of such solutions actually satisfy a linear 
differential equation of order five and six respectively provided the coefficients 
are sufficiently smooth. Indeed, it is the case that given any solution of a linear 
second order equation its m-th power satisfies an effectively computable linear 
differential equation of order m+1. 



1. The differential equations 

Consider the real linear second order differential equation in an unspecified in- 
terval /, of the form 

y" = p(x)y' + q(x)y, x € I, (1) 

where p, q are sufficiently smooth functions to allow the operations below to proceed. 
For example, if / is a solution of (TT]) and we are seeking the differential equation 
satisfied by the function y(x) — f(x) m , m > 3, then p,? e C m_1 (I). In what 
follows two linearly independent solutions of ([1]) will always be denoted by /, g. 
The case m = 2 is well-known and admits a myriad of applications: The result can 
be stated simply as saying that the quantities f 2 , fg,g 2 form a basis for a related 
linear third order differential equation, that is, cf., [[5], p. 298] 

y'" = 3p y" + (4q - 2p 2 + p')y' + 2(q'- 2pq) y, x £ I, (2) 

Its early appearance can be found in the books by Ince [4] and Whittaker and 
Watson [5], although the latter credit Appell [I], the case m = 2 appears to be 
due to Clausen [2J, see Gasper and Rahman [[3], p.xv]. In the developments stated 
below all proofs are lengthy but straightforward and can be ascertained by direct 
verification. The third power analog of this result is as follows: Whenever f'(x) ^ 
and g'(x) ^ 0, we have the following results: 
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Theorem 1.1. The functions f 3 , f 2 g, fg 2 , g 3 form a basis of the linear fourth order 
differential equation 

-6py"' + (lip 2 - V - I0q)y" + (7 PP ' + 30pq-6p 3 -p" - 10q')y' 
+ (9q 2 + 15pq' + 6 P 'q - 18p 2 q - 3g") y = 0. (3) 

The case m = 3 above is due to Appell [[T],p.213] but is not commonly known; 
none of the classic authors appear to mention it or any other generalizations. We 
have not been able to find any references for the next two cases, nevertheless they 
are interesting enough to be listed since the study of orthogonal polynomials has 
reached equations of degree at least six these days. 

Theorem 1.2. The functions / 4 , f 3 g, f 2 g 2 , fg 3 ,g A form a basis of the linear fifth 
order differential equation 

y {v) + c 4 y {lv) + c 3 y'" + c 2 y" + c iy ' + c o y = 0. (4) 

where 

C4 = — lOp, 

c 3 = 35p 2 - 10p' - 20q, 

c 2 = App' - 5p" - 30q' + 120p<? - 50p 3 , 

ci = 64q 2 + 24p 4 + h&p'q + Upp" - 46p 2 p' - 2Q8p 2 q + I28pq' - p'" + 7p' 2 - 18q" , 
c = 36pq" - I28pq 2 - 4q"' + 28p'q' - 104pY + Mqq' - 80pp'q + 8p"q + 96p 3 q. 

Theorem 1.3. The functions f 5 , f 4 g, f 3 g 2 , f 2 g 3 , fg 4 , g 5 form a basis of the linear 
sixth order differential equation 

y {m) + c 5 y^ + c 4 y {lv) + c 3 y'" + c 2 y" + c 1 y / + c y = 0. (5) 

where 

c 5 = -15p, 

c 4 = 85p 2 - 35q - 20p', 

c 3 = 350qp - 225p 3 - 70q' - 15p" + 165pp', 

c 2 = 274p 4 -63g"-1183p 2 g+52p /2 + 525pg'+81pp"-421pV + 259 g 2 + 266gp'-6p"', 

ci = 315p<f - 120p 5 - 28q"' - - I071pp'q + 326p 3 p' - 1295pg 2 + 1540p 3 g - 

lOlpV + 16pp"' + 25p'p" + 9lqp" + 266pY - 1183pY - 127pp' 2 + 518qq', 

c = U50q 2 p 2 -l295qq'p-310q 2 p'-5q ( ■ ^v '>+860p 2 p'q-535pp'q / -150pp"q+770p 3 q , - 

600p A q-100p' 2 q-355p 2 q"+45p"q'+80p'q"+10p'"q+155qq"+7Qpq'"+130q' 2 -225q 3 . 

Remark An (m+l)-th order linear ordinary differential equation with the basis 
elements f m , f m ~ 1 g, . . . , f l g^ , ■ ■ ■ , g m where i + j = m, can be derived as well. 
Thus, higher order equations of the type presented above can be integrated easily. 
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